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We study the equation of state of neutron matter using a family of unitarity potentials all of 
which are constructed to have infinite ^So scattering lengths Us- For such system, a quantity of 
much interest is the ratio ^ = Eq/Eq^'^'^ where Eq is the true ground-state energy of the system, 
and Eq^"^" is that for the non-interacting system. In the limit of Os —5- ±oo, often referred to as the 
unitary limit, this ratio is expected to approach a universal constant, namely ^ ~ 0.44(1). In the 
present work we calculate this ratio ^ using a family of hard-core square- well potentials whose as can 
be exactly obtained, thus enabling us to have many potentials of different ranges and strengths, all 
with infinite as. We have also calculated ^ using a unitarity CDBonn potential obtained by slightly 
scaling its meson parameters. The ratios ^ given by these different unitarity potentials are all close 
to each other and also remarkably close to 0.44, suggesting that the above ratio ^ is indifferent to 
the details of the underlying interactions as long as they have infinite scattering length. A sum-rule 
and scaling constraint for the renormalized low-momentum interaction in neutron matter at the 
unitary limit is studied. The importance of pairing in our all-order ring diagram and model-space 
HF calculations of neutron matter is discussed. 

PACS numbers: pacs 
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I. INTRODUCTION 

The 'unitary limit' of an ultracold Fermi gas refers to 
the special scenario where the inter-atomic interaction is 
tuned to have its scattering length as approaching infin- 
ity (os >> kp^), leaving the Fermi momentum kp &s the 
only relevant length scale in the many-body system. This 
limit was discussed as early as 1999 by Bertsch [l|, who 
challenged many-body theorists with the question "What 
are the ground state properties of neutron matter, inter- 
acting with an infinite scattering length?" . Under such 
condition the Fermions are strongly interacting, thus the 
determination of the ground state is highly non-trivial. 
However, the equation of state at this limit is expected to 
be an universal expression Eq = S^Eq^'^'^ where ^ is an uni- 
versal constant for any underlying Fermion systems. It is 
of much interest to determine ^ experimentally or derive 
it theoretically, many such attempts having been made 
P-[27j. By far, the best numeric estimate on is consid- 
ered to be C = 0.44(1) 0, 0.42(1) ^ and 0.40(1) [H 
according to quantum Monte Carlo methods. The scat- 
tering length of trapped ultracold Fermi gases can in fact 
be magnetically tuned [2§| by way of the Feshbach reso- 
nance. This is indeed a very important achievement, and 
it is based on this that the above constant ^ for trapped 
ultracold Fermi gases can be experimentally measured. 
Several recently reported experimental values for ^ are 
0.46 ±0.05 0, OAGto ll S and 0.39(2) to 0.435(15) 0. 

In a previous work j2l|, we have calculated the con- 
stant ^ for neutron matter using a tuned CDBonn re- 
alistic nucleon-nuclcon (NN) potential [30]. Unlike the 
case of trapped Fermion gases, here the scattering length 
of NN interactions still can not be tuned experimentally. 



Thus we have chosen to tune the NN interaction by ad- 
justing its meson parameters, in line with the Brown-Rho 
scaling [32| that the meson masses in nuclear medium is 
suppressed compared with their in-vacuum values. The 
^5*0 scattering length of the original CDBonn potential is 
as=-18.9 fm, which is already fairly long. We have found 
that to have a very large scattering length (such as -12070 
fm) the meson parameters of the CDBonn potential need 
to be adjusted only slightly (about 2%) [2_l|. This tuned 
CDBonn potential actually has given ^ quite close to 0.44 
for a wide range of densities (~ 0.02 to ^ 0.09 fm~^) [2l|. 

In the present work, we would like to investigate mainly 
the following question: Is the constant ^ given by other 
unitarity potentials, defined as those with infinite scatter- 
ing lengths, also close to 0.44? This constant is supposed 
to be an universal constant, and then its value should 
be independent of the detail structure of the potentials 
as long as they have infinite scattering length. In other 
words, ^ given by all other unitarity potentials should in 
principle be the same. It may be difficult to prove this 
analytically. Before one can do so, it should be useful 
and of interest to first check this universality property 
numerically as we shall do in the present work. As it is 
rather involved to tune the various realistic NN poten- 
tials [sil - lssj to infinite scattering length, we have chosen 
to calculate ^ from a family of simple hard-core square- 
well (HCSW) potentials 



V{r) = Vc, r < rc 

Vb] rc <r <rb 
0; r > Tfc. 



(1) 



An advantage of this type of potentials is that their scat- 
tering lengths can be given analytically and thus exact 
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unitarity potentials can be readily obtained. The above 
HCSW potentials are clearly very different from the re- 
alistic NN potentials, yet as we shall discuss later the 
constant given by various HCSW potentials are all amaz- 
ingly close to 0.44, as is also given by the tuned CDBonn 
potential. 

In the following we shall present some details of our 
calculations. In section IIA, we shall first briefly outline 
the renormalized low-momentum interaction Viow-k (ssl - 
l40j on which our calculations will be based. We shall 
carry out calculations using both HCSW and CDBonn 
potentials, both having strong repulsive cores. As is well- 
known, these hard-core potentials are not suitable for 
being used directly in many-body calculations. In con- 
trast, the low-momentum interaction Viow~k is a smooth 
and energy-independent potential, and it has been ex- 
tensively applied to calculations of nuclear matter [4l|, 
nuclear structure [i^l and neutron star i45]. This inter- 
action is obtained using a renormalization group (RG) 
method where the high-momentum components of an un- 
derlying potential beyond a decimation scale A are inte- 
grated out. The low-energy phase shifts are preserved by 
the RG procedure and so is the scattering length. In sec- 
tion IIB the procedures for constructing our model uni- 
tarity potentials are presented. Analytical expressions 
for the scattering length of HCSW potentials are given. 
From these expressions, one can construct an unlimited 
number of unitarity potentials (which all have infinite 
scattering length). 

In section IIC, we shall describe how we calculate 
the neutron matter equation of state (EOS) from the 
low-momentum interaction Viow-k- Two methods will 
be employed: a ring-diagram method [2l|, |4l| and a 
model-space Hartree-Fock (MSHF) method. There have 
been several many-body methods for calculating the 
EOS for large (number of particles approaching infinity) 
Fermionic systems such as atomic gases and nuclear mat- 
ter. For a long time, the Brueckner-Hartree-Fock (BHF) 
theory [46!-i48i] was the commonly used framework for 
such EOS calculations; an inconvenience of this method 
is that the BHF G-matrix interaction is energy depen- 
dent, adding numerical complications. We shall use the 
energy independent Viow-k interaction mentioned earlier 
for the EOS calculations, to circumvent the above incon- 
venience. With the use of the V5ou,_fc interaction, the 
MSHF method provides a simple way for calculating the 
constant £,. In fact MSHF is rather similar to BHF: In 
BHF the particle- particle correlations above kp are incor- 
porated by way of the energy-dependent G-matrix, while 
in MBHF they are included via the energy-independent 
Viow~k interaction. As we shall discuss later, the con- 
stant ^ given by the ring-diagram method and the MSHF 
method are in fact nearly identical (both close to 0.44). 
We shall also discuss that the MSHF mean field poten- 
tial should obey certain special constraint at the unitary 
limit. The importance of pairing in fermionic systems 
such as neutron matter has been well known and exten- 
sively studied (see e.g. [13, H^l and references quoted 



therin). In section IIC, we shall discuss that pairing has 
also played an important role in our all-order ring dia- 
gram calculations where a pair of fermions are allowed 
to interact any number of times, forming a quasi-boson 
composed of a coherent pair of fermions. Our results will 
be presented and discussed in section HI. A summary 
and conclusion will be given in section IV. 



II. FORMALISM 

A. Low-momentum interactions 

Since our calculations are largely dependent on the re- 
cently developed low-momentum interaction Viow-k, it 
may be useful to present first a short review on this in- 
teraction. As we know, the strong short-range repulsion 
contained in Vnn usually requires a special renormaliza- 
tion treatment before being used in many-body calcula- 
tions. A familiar such treatment is the BHF G-matrix 
method where the G-matrix interaction is obtained by 
summin g th e particle-particle ladder diagrams to all or- 
ders |46l - t48| . An inconvenient aspect of this method is the 
energy dependence of G, making the many-body calcula- 
tions based on G rather complicated. In the past several 
years, there has been much progress in the RG approach 
to the NN interaction. A central idea here is that to de- 
scribe the low-energy properties of a physical system it 
should be adequate to employ only a low-momentum ef- 
fective interaction confined within a momentum decima- 
tion scale A. Viow~k is such a low-momentum effective 
interaction. Starting from a realistic NN potential VnNi 
Viow-k is obtained from the following T-matrix equiva- 
lence equations [35l - l40| : The half-on-shell T-matrix for 
the full-space interaction is 



r(/c',fc,fc2) = VNN{k',k) 



TT Jo 



^2 _ 2 9 (2) 



where V denotes principal-value integration and the in- 
termediate state momentum q is integrated from to 
oo covering the whole space. Then we define a low- 
momentum half-on-shell T-matrix by 



Tiow—k{k ,k,k ) — Viow—k{k , k) 



n 



^ Viow-k{k', q)Tio„_k{q, k, fc^) 

k^-q' 



q^dq, (3) 



where the intermediate state momentum is integrated 
from to A, the model space cut-off. The low momen- 
tum Viow-k is then obtained by requiring the equivalence 
of T-matrix in the model space, 

T(fc', fc, fc2) = Tio^,-k{k\ k, fc2); (fc', k) < A. (4) 

The low-momentum interaction Viow-k is obtained by 
solving the above three equations. Since this proce- 
dure preserves the half-on-shell T-matrix, it certainly 
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preserves the low-energy (< A^) phase shifts and scat- 
tering length Us of V/vat- This ensures that if Vatjv is 
tuned to be a unitarity potential (a^ — )■ ±00), so is the 
corresponding Viow-k- 

B. Hard-core square- well unitarity potentials 

There are several realistic NN potentials (CDBonn[30|. 
Argonne[3ll|, Nijniegen[32|. Chiralfss'l) which all describe 
the two-nucleon low-energy experimental data very accu- 
rately. It would be of interest to obtain unitarity poten- 
tials from them by tunning their interaction parameters 
slightly. But technically this is not easy to carry out. So 
far only the CDBonn potential has been tuned to attain 
this limit (a3=-12070 fm)j21,]. In the present work, we 
choose hard-core square- well (HCSW) potentials as given 
by Eq.(l) for further studying the properties of neutron 
matter at unitary limit. An advantage of using them is 
that their scattering length can be analytically ob- 
tained, allowing us to study the properties of neutron 
matter using many HCSW potentials with any chosen 
scattering lengths (including infinity). 

As indicated in Eq.(l), the HCSW potential is charac- 
terized by the height Vc of the repulsive core, the depth 
Vb of its attractive well, and the respective ranges Tc and 
Tf,. In the present work, we consider neutron matter with 
interactions only in the ^Sq channel, whose phase shift S 
is readily obtained from Eq.(5) and (6). Namely 

K3 

tan{6 + K^rb) = —tan{K2rb + a) (5) 
K2 

tan{a + K2rc) = -—tanh{Kirc) (6) 
Ki 

with 

/ JTl 

A-2 = ^/(E-Vi)^, 

where E is the scattering energy in the center-of-mass 
frame. From the above results, the scattering length 
is obtained from a low energy expansion of kcotS as 



with 

A = K10K2Q - K2Qtanh{KiQrc)tan[K2oirb - Tc)], 
B = K2otanh{Kiorc) + Kwtan[K2o{rb - r^)] - rbKioK2o 
+rbK2Qtanh{Kiorc)tan[K2o{rb -^c)], (8) 

where 



TABLE I: Three different unitarity HCSW potentials. 



Potentials K[MeV] rc[fm] 


V6[MeV] 


rt [fm] 


a^xlCfm] 


re [fm] 


HCSWOl 3000 0.15 


-20 


2.31 


15.2 


2.36 


HCSW02 3000 0.30 


-30 


2.03 


3.38 


2.21 


HCSW03 3000 0.50 


-50 


1.81 


-4.58 


2.20 



The effective range Tg for this potential can also be de- 
rived analytically; the result is not presented here be- 
cause it is fairly lengthy. 

Eq.(7) implies that the condition for being a unitarity 
potential (infinite scattering length) is j4=0, namely 

rb — Vc = tan~^[ ^ -]. (9) 

K20 K2atanh{Kiorc) 

It is well known that when the potential is tunned to the 
unitary limit (the Feshbach resonance), it has a bound 
state with its energy approaching to zero. It is readily 
checked that the condition for having such a bound state 
is the same as the ^=0 one given above. 

The above condition enables us to construct any num- 
ber of unitarity potentials by varying any three of the four 
parameters Vb, Tc and r^. Some sample such poten- 
tials are listed in Table I. As seen, they all have very large 
(infinite) scattering lengths while the potentials them- 
selves are significantly different from each other. For ex- 
ample, the depth of the attractive parts of them changes 
from -20 to -50 MeV. The effective ranges of these po- 
tentials are also listed, all being close to the unitarity- 
CDBonn's r,, of 2.54/m~^. These potentials will be used 
to calculate the universal ratio ^ = Eq/Eq^'^'^ for neutron 
matter. As we shall soon discuss, the ratios ^ given by 
these three largely different potentials are in fact nearly 
identical to each other. 



C. Ring-diagram and model-space Hartree-Fock 
methods 



0--0 




(a) (b) (c) 

FIG. 1: Diagrams included in the all-order pphh ring-diagram 
summation for the ground state energy shift of nuclear matter. 
Each dashed line represents a Viom-k vertex. 

We use the ring-diagram method [2l|, |4l| to calculate 
the neutron matter EOS. In this method, the ground- 
state energy shift AE is given by the all-order sum of 
the pphh ring diagrams as illustrated in Fig.l, where (a), 
(b) and (c) are respectively 1st-, 4th- and 8th-order such 
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diagrams. (A£;o is defined as (Eq - E},"'"'') where Eq 
is the true ground state energy and E^^'^^ that for the 
non- interacting system.) In this way, we have 

•^0 m ijkKK 

xy,;(fc/,A)(zj|v^o,^_fc|fcO, (10) 

where the transition amplitudes are 

Y:^{kl,\) = (*,„(A,A-2)|a/afe|«'o(A,A)), (11) 

where ^o(A, A) denotes the true ground state of nuclear 
matter which has A particles while 4'm(A, A — 2) the mth 
true eigenstate of the (^-2) system. Note that A is a 
strength parameter, integrated from to 1. Using the 
above Ai?o, the ratio is readily obtained, namely 1 
+l^E^lEl''\ 

By summing the pphh ring diagrams to all orders, the 
above amplitudes Y are given by the coupled equation of 
the form 

AX + BY = LoX 
A*Y + B*X = -ujY, (12) 

where X is y's RPA partner. [i^l w is the excitation en- 
ergy corresponding to either {Em{A — 2) — Eq{A)) or 
{Em{A + 2) — Eo{A)). It is readily recognized that the 
above RPA equation is of the same form as the well- 
known quasi-boson RPA equations (isj . resulting from 
treating the many-body system approximately as a col- 
lection of quasi bosons (each composed of a pair of inter- 
acting fermions). Thus our all-order ring-diagram calcu- 
lation is actually a quasi-boson approximation for neu- 
tron matter, treating it as a system of quasi bosons. 
(This is not unexpected as RPA is a quasi-boson ap- 
proximation.) Our all-order ring-diagram calculation is 
practically a pairing calculation in the sense that we in- 
clude only the boson-forming ring diagrams to all orders 
(both forward and backward ladder diagrams for a pair of 
fermions interacting any number of times) , allowing them 
to form a coherent pair (quasi-boson). Note, however, 
that our calculation becomes HF if we include only the 
Ist-order ring diagram, because in this case the above Y 
amplitude becomes Y[ij) = riirij where ni=l for i < kp 
and =0 otherwise. Since we include only the ^Sq in- 
teraction, the quasi bosons in our present ring-diagram 
calculation are in fact all BCS-pairing bosons (^5*0). 

The ring-diagram method described above is a model- 
space formalism, where all nucleons are confined within 
a momentum model space P(fc < A), and the decimation 
scale for the Viow-k interaction is the same A. The ex- 
pression Eq.(lO) for the energy shift AEq is rather com- 
plicated, and is not convenient for studying the proper- 
ties of the underlying interaction at the unitary limit. 
To have a clearer way to study these properties, we have 
considered a simpler method to calculate AEq, namely 



a model-space HF method (MSHF). In this method, we 
have 

AEo = l^k,M<kAkik2\V,t-k\kik2). (13) 

This is also a model-space approach where all nucleons 
are confined within a model space P{k < kp), but MSHF 
has a special feature that the interaction used is renor- 
malized according to the same P{k < kp), namely the 
Viow-k interaction is calculated with A = kp as denoted 
by Vi^^_^. A main difference between MSHF and the 
usual HF calculations is that the interaction used in the 
latter is not required to be renormalized according to the 
above model space while it is so for the former. In fact 
MBHF is equivalent to the ring-diagram method when 
one chooses A to be its smallest value allowed by the 
Pauli principle, namely kp. When A — kp, diagrams 
like (b) and (c) of Fig.l no longer exist (since particles 
with k > kp are disallowed in the model space); only 
diagram (a) remains, which is just MBHF. It may be 
mentioned that MBHF is rather similar to the famil- 
iar BHF method [46l - l48l |. They both employ the same 
model space P{k < kp) (all nucleons confined within 
the Fermi sea), but they employ different model-space 
effective interactions. In BHF the energy-dependent G- 
matrix interaction is used while in MBHF the energy- 
independent Viow-k interaction is employed; both include 
renormalizations from all-order particle-particle correla- 
tions beyond kp, but for G the renormalization is carried 
out using an energy-dependent method while an energy- 
independent one is applied for Viow-k- The MBHF 
ground-state energy is given by the sum of the kinetic 
energy and a simple integral, namely 

Eo 3 8 f^" ,2^, 3fc fc3 

T^r^+^X ^ ^^[1 - 2fc; + 24] 

xY^{2J^^\){a,k\Vt,_^\a,k). (14) 
n^k^ 

In the above £p= „ , k denotes the two-nucleon rela- 
five momentum and Jq. is the total angular momentum 
for the partial wave a. (We consider neutrons with inter- 
actions only in the ^Sq partial wave, and thus the above 
summation has only one term.) The above result may be 
useful for studying the properties of the Viow-k interac- 
tion at the unitary limit. If the ratio Eq/Eq^'^'^ is equal 
to a constant ^ at this limit, then the above equation 
implies that at this limit must satisfy 

xeSo,k\V,';-_,\'So,k). (15) 

This is a rather strong constraint for the interaction. 
Clearly there are many potentials which can satisfy the 
above constraint, allowing many potentials to have the 
same ratio ^. In the following section, we shall discuss 
and numerically check this constraint. 
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Within the above MSHF framework, the single-particle 
(s.p.) potential U is given [42| as 

1 fi Z"*^" 

Uih) = ^(2J„ + 1){- / edk{ka\Vit-k\ka) 
„ ^ Jo 



III. 



RESULTS AND DISCUSSIONS 



/ ^ kdk[kj, ~kl + 4fc(fci - fc)] 

with 

fc_ = (fci^ - fci)/2, fc+ = (fcF + fci)/2. 
The MSHF s.p. spectrum is 

h^kl 



(16) 



e(ki 



2m 



U{k, 



which can be well approximated by 



(17) 



(18) 



where to* is the effective neutron mass in medium and 
A is a constant related to the depth of the potential 
well. The MSHF ground-state energy can be expressed 
in terms of m* and A. Then if this energy is equal to 
^eI''"" at the unitary limit, the MSHF m* and A should 
satisfy the linear constraint 



1 TO 5A 
^=2 + 2^^+6^- 



(19) 



In the following section, we shall also discuss and numer- 
ically check the above linear constraint concerning the 
MSHF sp potential. 
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FIG. 2: Ring-diagram unitary ratios ^ given by different uni- 
tarity potentials. 



We first calculate the ratio ^ with the ring-diagram 
method (see Eq.(9)) using both the three HCSW unitar- 
ity potentials of Table I and the unitarity CDBonn poten- 
tial (a3=-12070 fm) [IJ. As shown in Fig.2, the results 
given by these four potentials are all quite close to each 
other and to 0.44. A common decimation scale of A = 2.1 
fm^^ has been used for the results presented there. We 
have found that the f values given by the ring-diagram 
calculations using the HCSW potentials are rather sta- 
ble (variation less than ^ 0.002) with A in the range be- 
tween ~ 2.0 and ^ 2.4 fm~^. As found in ^21'], they are 
also quite stable for the unitarity CDBonn potential in 
the same range. Since these potentials have given nearly 
identical results for the ration ^, it is of interest to com- 
pare their V^Nik, k) matrix elements from which the ^s 
are calculated. Such a comparison is presented in Fig. 3, 
and as seen they are actually quite different. It is indeed 
surprising that the ratios ^ given by these vastly differ- 
ent interactions are nearly the same. To illustrate the key 
role of the unitary limit, we further calculate the ratios ^ 
near the unitary limit as displayed in Fig. 4. When Ug is 
away from the unitary limit, the f values given by these 
potentials are noticeably different, but they converge to 
a common value only as l/flg approaching zero. Above 
results strongly suggest that the ratio ^ is indifferent to 
the details of the underlying potentials, as long as they 
have infinite (or very large) scattering lengths, namely 
they are all unitarity potentials. 




CDBonn 
HCSW01 
HCSW02 
HCSW03 



k[fm"^] 



FIG. 3: Diagonal matrix elements VNN{k,k) of different uni- 
tarity potentials. 

Using the above four unitarity potentials, we have also 
calculated the ratios ^ using the MSHF method (see 
Eq.(14)). The results are presented in Fig. 5, and as seen 
the results given by the four potentials are nearly identi- 
cal and they are all quite close to 0.44 over a wide range of 
densities. It is of interest to note that the MSHF results 
shown are remarkably close to the ring-diagram ones of 
Fig.2. This close agreement is an indication that the hole- 
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0.2 0.4 0.6 0.8 



k[fm"^] 



FIG. 4: The ring-diagram unitary ratios ^ near the unitary 
limit, kp = 1.2 fm~^ is used. 



FIG. 6: Diagonal matrix elements of Vjj,^_j, from different 
unitarity potentials. A = fc^ = 1-2 fm~^ is used. 



hole correlation diagrams like diagram (c) of Fig.l may 
not be important for pure neutron matter. The contribu- 
tion from the particle-particle correlations like diagrams 
(a) and (b) is included in both ring diagram and MSHF 
calculations, but diagrams like (c) are included only in 
the former. It may be noted that in both Fig. 2 and Fig. 5 
our results are nearly constant for kp between ^0.9 and 
1.4 fm"""^ (corresponding to density range ^ 0.016 to 
^ 0.092 fm^''). Outside this range, small variations of ^ 
start to appear. These variations are an indication that 
the methods we have employed for the calculation is not 
adequate for the high and very low density regions in- 
dicated above. So far we have employed only two-body 
Viow-k interactions in the calculation. Three body in- 
teractions have not been considered in the present work, 
and their effects may be important. Further studies in 
this direction are needed. 
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FIG. 5: MSHF unitary ratios ^ given by different unitarity 
potentials. 

In the MSHF method, the ratio ^ is calculated from 
the model-space effective interaction VJo^_j,. We have 
seen from Fig. 5 that the ratios ^ given by the four dif- 



ferent potentials are nearly the same. Does this mean 
that their are also nearly the same? In Fig. 6 the 

matrix elements of this interaction for four different uni- 
tarity potentials are presented. As seen, they are clearly 
not identical to each other; there are significant differ- 
ences among them. Eq.(14) is in fact a strong constraint 
for the above interactions. To have ^ being a univer- 
sal constant, the left-hand side of Eq.(14) is proportional 
to hp. Consequently the integral in the right-hand side 
must be equal to C x fc^, C being a constant, noting that 
its integrand and integration limit are both dependent on 
kp. This implies that at the unitary limit the V^o^_j, in- 
teraction must satisfy certain stringent requirements so 
as to satisfy the above constraint. 

Similar to the effective interactions commonly used in 
effective field theories [29, Esl, we have found that our 



1^ 



interactions can be highly accurately simulated 



by low-order momentum expansions of the form 



ki 



'fee 



(20) 



where Vq, V2 and V4 are constants (independent of k but 
dependent on kp). (The rms deviations for fitting the 
results of Fig. 6 by the above expansion are ~ 2 x 10~^ 
and ^ 2 X lO"-^ respectively for the CDBonn and HCSW 
cases, the fitting being very good.) In terms of these 
constants, Eq.(14) assume a rather simple form, namely 



Y^(e-l) = fcF(y 



10 



3^4 



(21) 



The above is an interesting sum-rule and scaling con- 
straint, namely at the unitary limit the strength sum 
(Vo/3-l- V2/2-I-3V4/7O) is a constant for any kp (density) 
and, in addition, it scales with 1/kp (i.e. proportional 
to 1/kp). In Table II we present some sample results, to 
check how well do they satisfy the above constraint. As 
seen our results satisfy this constraint very well. For each 
kp, the values of the above sum given by the four differ- 
ent unitarity potentials are all quite close to each other. 
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The values of this sum for different kp are different but 
they are all giving f close to 0.44, in close agreement with 
the 1/kp scaling. 

TABLE 11: Low-order momentum expansion 

of^;o»-fc- Listed 

are the coefficients Vb, V2 and V4 of Eq.(20), with the sum 
(Vo/3 + 1/2/10 + 314/70) denoted as Sum. Four unitarity 
potentials are used. 



Vivjv kp 


[fm-^] 


Vi,[fm] 


1^2 [fm] 


14 [fm] 


Sum 




CDBonn 


1.2 


-2.053 


3.169 


-1.801 


-0.445 


0.434 


HCSWOl 




-2.001 


2.865 


-1.402 


-0.441 


0.439 


HCSW02 




-1.904 


2.373 


-0.999 


-0.440 


0.440 


HCSW03 




-1.893 


2.261 


-0.825 


-0.440 


0.439 


HCSWOl 


1.0 


-2.102 


2.202 


-1.070 


-0.526 


0.442 


HCSWOl 


1.4 


-1.945 


3.584 


-1.983 


-0.375 


0.443 



-0.68 




-0.88 I — ■ — ' — ■ — ' — ■ — ' — ■ — ' — ■ — ' — ■ — I 
1.05 1.1 1.15 1.2 1.25 1.3 1.35 

m/m* 

FIG. 7: MSHF effective mass m* and well depth A at the 
unitary hmit. Results calculated from the unitarity potentials 
are denoted by 'star', 'dot' and 'triangle' respectively for kp 
in the ranges (0.8-0.9), (0.9-1.4) and (1.4-1.5) fm'^ The solid 
line represents the linear expression Eq.(18) with ^=0.44. See 
text for other explanations. 

By far, the best numeric estimates of the ratio ^ 
have been obtained from the variation al q uantum- Monte- 
Carlo (QMC) calculations (see e.g. [2J| and references 
quoted therein). It would be useful and of interest to 
compare our ring-diagram and MSHF calculations with 
the QMC ones. As reviewed in [131, superfluid-QMC 
calculations are generally based on the variation princi- 
ple (5[($BCs|^'ji?*j|*BCs)/($BCs|*j*j|^BCs)] = 0, 

where ^bcs is a BCS-paired trial wave function of defi- 
nite number of particles and is a r-space Jastrow cor- 
relator whose inclusion depends on the interaction em- 
ployed; its inclusion is necessary when an interaction 
with a 'hard-core' repulsion is employed |23| but not 
so for smooth potentials without 'hard core^a. flTI. [2^ . 
Clearly pairing is incorporated in superfluid-QMC by us- 
ing a BCS-paired trial wave function. As discussed in 
section II, we emphasize pairing differently, namely by 
including only the boson-forming ring diagrams to all or- 
ders (both forward and backward ladder diagrams for a 



pair of fermions interacting any number of times), and 
in this way our all-order ring-diagram calculation treats 
the system as a collection of quasi bosons. Despite this 
difference, it is encouraging that our result of ^ ~0.44 is 
satisfactorily close to the superfluid-QMC value of ~ 0.40 
to ^ 0.44 [2^ [23|. It may be pointed out that when 
we include only the first order ring diagram, our ring- 
diagram calculation reduces to the usual HF calculation. 
For A — 2.1/m~^ and kp = 0.8/m~^ we have carried 
out HF calculations (i.e. including only the Ist-order 
ring diagram) and obtained £,hf —0.55, for all the four 
unitarity potentilas mentioned earlier. This S^hf is sig- 
nificantly larger than the corresponding all-order-ring re- 
sult of 0.44. Normal-QMC calculations (where the trial 
wave function is a closed Fermi sea) have given ^ ~0.54 
[l^illl]- We believe that our Ist-order ring calculation 
is similar to that of the normal-QMC while our all-order 
ring-diagram one is similar to that of the superfluid-QMC 
(where the trial wave function is ^bcs)- Further studies 
of the above possible connections are needed. 

We have discussed in section II that at the unitary limit 
the MSHF s.p. potential has a special property, namely 
its effective mass m* and well depth A satisfy a linear 
constraint dependent on the unitary ratio ^ (see Eq.(18)). 
Here let us check how well is this constraint obeyed by our 
results. Similar to Fig. 5, we have calculated m/m* and 
A/ep using the CDBonn and HCSW unitarity potentials 
using a wide range of kp from 0.8 to 1.5 fm~^ . Our results 
are plotted in Fig. 7, and as seen most of our data are 
located near the straight line corresponding to ^=0.44. 
In the figure, the results for kp in the range of (0.9- 
1.4)fm~^ are quite close to the ^=0.44 straight line. For 
kp outside this range, they are slightly off the line, which 
is consistent with the results of Fig. 5. 

It is promising that at the unitary limit the above 
MSHF framework has given rather satisfactory results 
for the ratio ^ (^ 0.44), in close agreement with both the 
superfluid-QMC results [H, H^l and those from our all- 
order ring diagram calculations. In addition, this method 
is fairly simple and transparent; as indicated by Eq.(15) 
^ is given by a simple integral and the HF potentials 
at the unitary limit satisfy the constraint of Eq.(19). It 
may be mentioned, however, that this MSHF is only a 
limited effective theory; it is the A = kp limit of the 
ring-diagram formalism. In this limit the model space 
is one dimensional, and the model-space ground state 
wave function is \kp), the closed Fermi sea. The BCS 
pairing gap Abcs is an important quantity, but this in- 
formation is not contained in the MSHF effective theory. 
Abcs is the lower bound for the lowest excitation en- 
ergy of the many-body system. MSHF, however, is a 
one-dimensional effective theory which can describe only 
the ground state of the system; it provides no information 
about excited states. The effective MSHF theory is not 
capable to provide information about AbcSj even though 
the MSHF wave function \kp) is completely {k t, —k l) 
paired. To calculate Abcs^ one may need to employ 
the number-nonconserving Green's function framework 
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[sol ISlj which has been commonly used in supercoduc- 
tivity (superfluidity) calculations. (The Green's function 
method used in the present work is based on a number- 
conserving formalism, and is not suitable for calculating 
^BCS-) Using several Skyrme effective interactions, Su. 
et al. (sil l have performed such number-nonconserving 
BCS calculations for nuclear matter. We plan to repeat 
their BCS calculations using instead the RG Viow-k in- 
teractions, and will report our results in a separate pub- 
lication. 



IV. SUMMARY AND CONCLUSION 

Using several different unitarity potentials, defined as 
having infinite (very large) scattering length a^, we have 
calculated the ground-state energy ratio ^ = Eq/Eq^'^^ 
for neutron matter over a wide range of densities. A 
main purpose of our study was to check if the ratios so 
obtained are 'universal' in the sense that they have a 
common value (or nearly so), independent of the details 
of the potentials as long as they have as — t- ±oo. We 
have used four unitarity potentials: One is the unitarity 
CDBonn potential whose meson exchange parameters are 
slightly tuned so that its becomes very large (-12070 
fm). The other three are square- well 'box' potentials with 
both hard-core repulsion and exterior attraction. The 
Us of these box potentials can be obtained analytically, 
making it easy to construct many unitarity box potentials 
of different ranges and depths. 

We have calculated ^ using two methods: a 
ring-diagram method and a model-space Hartree-Fock 
(MSHF) method. An important step in both methods is 
the transformation of the unitarity potentials, which have 
strong short-range repulsions, into low-momentum inter- 
actions Viow-k which are smooth potentials, convenient 
for many-body calculations. The transformation is car- 
ried out using a renormalization group method which pre- 
serves low-energy phase shifts and the scattering length. 



We have found that the ratios f given by the ring-diagram 
and MSHF methods are practically identical over a wide 
range of densities (from 0.02 to ~ 0.09 fm^"^). This is 
an indication that the particle-particle hole-hole correla- 
tions (as shown by diagram (c) of Fig.l) are not impor- 
tant for neutron matter. Such correlations are included 
in the ring-diagram method but not in MSHF. The effect 
of the particle-particle ladder correlations (as shown by 
diagram (b) of Fig.l) are included in both. The MSHF 
method is considerably simpler than the ring-diagram 
one. It may provide a promising method for neutron 
matter. Further study of this method for neutron matter 
may be useful and of interest. 

The CDBonn and the three box unitarity potentials are 
very different from each other, but the ratios ^ predicted 
by them using either of the above methods are all remark- 
ably close to 0.44, over a wide range of densities. This is 
indeed a rather interesting and surprising result, indicat- 
ing that ^ may be perceived a universal constant. For the 
MSHF case, this result is related to some special prop- 
erties of the Viow-k interaction at the unitary limit. We 
have found that the Viow-k interaction for MSHF can be 
accurately simulated by a low-order momentum expan- 
sion of the form [Vq + V2{k/ kp)'^ + V4{k/ kp)'^], and at the 
unitary limit the strength sum (Vo/3 -I- V2/IO -I- 3V4/7O) 
calculated from the above potentials satisfies the con- 
straint of Eq.(20) very well, for a wide range of kp. It 
is also found that at the unitary limit the parameters 
TO* and A of the MSHF mean field potentials given by 
the above potentials all obey a linear constraint satisfac- 
torily. In conclusion, we believe that our results provide 
strong numerical support to the conjecture that the ratio 
^ is universal. 
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